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ABSTRACT 



Aims. Our aim is to understand the evolutionary sequence of high-mass star formation from the earliest evolutionary stage of high- 
mass starless cores, via high-mass cores with embedded low- to intermediate-mass objects, to finally high-mass protostellar objects. 
Methods. Herschel far-infrared PACS and SPIRE observations are combined with existing data at longer and shorter wavelengths to 
characterize the spectral and physical evolution of massive star-forming regions. 

Results. The new Herschel images spectacularly show the evolution of the youngest and cold high-mass star-forming regions from 
mid-infrared shadows on the Wien-side of the spectral energy distribution (SED), via structures almost lost in the background emis- 
sion around lOOyt/m, to strong emission sources at the Rayleigh- Jeans tail. Fits of the SEDs for four exemplary regions covering 
evolutionary stages from high-mass starless cores to high-mass protostellar objects reveal that the youngest regions can be fitted by 
single-component black-bodies with temperatures on the order of 17 K. More evolved regions show mid-infrared excess emission from 
an additional warmer component, which however barely contributes to the total luminosities for the youngest regions. Exceptionally 
low values of the ratio between bolometric and submm luminosity additionally support the youth of the infrared-dark sources. 
Conclusions. The Herschel observations reveal the spectral and physical properties of young high-mass star-forming regions in detail. 
The data clearly outline the evolutionary sequence in the images and SEDs. Future work on larger samples as well as incorporating 
full radiative transfer calculations will characterize the physical nature at the onset of massive star formation in even more depth. 
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1. Introduction 

Characterizing the early evolutionary stages of (high-mass) star 
formation is difficult because the cold cores (-10 to 20 K) 
have the peak of the spectral energy distribution (SED) at far- 
infrared wavelengths which were hardly accessible until re- 
cently. Herschel has changed this situation completely, and we 
are now able to study the young (massive) star-forming regions 
in detail. For this early Herschel study we selected the com- 
plex associated with the High-Mass Protostellar Object (HMPO) 
IRAS 18223-1243 and the In frared Dark Cloud I RDC 18223-3 
at a distance of -3.7 kpc ( Sridh aran et all 120021) . This region 
hosts most evolutionary stages in massive star formation, starting 
with high-mass starless cores, continuing with high-mass cores 
with embedded low- to intermediate-mass protostars potentially 
forming massive stars, to HMPOs with already embedded and 
likely still accreting massive protostars (Fig.[T]). 

Previou s studies focused either on the HMPO IRA S 18223- 
1243 (e.g.. ISridharan et al.1 l2002t iBeuther et ail l2002h or at a 
high-mass core with embedded low- to intermediate-mass pro- 
tostar IRDC 18223-3 (JGarav et all l200J lBeuther & SteinackeiL 



120071: iFallscheer et alll2009l) . However, the whole complex with 
the youngest starless cores has not been studied so far. Here we 
combin e Herschel PACS and SPIRE data at far-infrared wave- 
lengths dPilbratt et all l2010HPoglitsch et alll20l5lGriffin etafl 
|2010|) with Spitzer mid-infrared data and (sub)mm continuum 
observations at the Rayleigh- Jeans tail of the spectrum. 



* Herschel is an ESA space observatory with science instruments 
provided by Principal Investigator consortia. It is open for proposals 
for observing time from the worldwide astronomical community. 



2. Data and observations 

The cl oud complex with a si ze of - 4' x 9 ' was observed with 
PACS dPoglitsch et all l2010h on Herschel dPilbratt et all l2010h 
on 2009 October 9 within the science demonstration program. 
Scan maps in two orthogonal directions with scan leg lengths 
of 18' and 6', respectively, were obtained with the medium scan 
speed of 20"/s. The raw data hav e been reduced with t he HIPE 
software, version 3.0, build 455 (JPoglitsch et all I201G) . Beside 
the standard steps leading to level- 1 calibrated data, a second- 
level deglitching as well as a correction for offsets in the detector 
sub-matrices were performed. Finally, the data were highpass- 
filtered, using a median window of the size of the full scan legs, 
to remove the effects of bolometer sensitivity drifts and the 1/f 
noise along the course of the data acquisition. We masked out 
emission structures (visible in a first iteration) before computing 
and subtracting this running median. Thereby over- subtraction 
of source emission in the highpass filtering step can be mini- 
mized. The flux correction factors provided by the PACS ICC 
team were applied. The beam sizes of the 70, 100 and 160 /mi 
data are -5.6", -6.8" and -11.4", respectively. The 100 yum 
PACS flux within an aperture of 40" toward the IRAS source 
18223-1243 is fully consistent with the calibration uncertainties 
of both missions (-17% measured difference). 

Maps at 2 50, 350, and 500 ^m were obtained with SPIRE 
dGriffinetalll2010h on 2009 October 19. Two 14' scan legs were 
used to cover the source. The data were processed within HIPE 
with the standard photometer script up to level 1 . During base- 
line removal, we masked out the high-emission area associated 
with the IRAS source. Because no cross-scan data were obtained 
for these observations, the iterative de-striping algorithm was in- 
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Fig. 1. The color-scale presents the wavelengths as marked in each panel. 70, 100 and 160 /mi are from PACS, 250/mi is from 
SPIRE. The contours show the SCUBA 850/mi continuum emission from 0.05 to 0.35 Jy beam -1 (in 0.1 Jy beam -1 steps) and 
additionally two stronger contours at 0.65 and 1.25 Jy beam -1 . The right panel labels which part of the region are dominated by 
what evolutionary stage. The circles and labels in the 160 /mi panel mark the sources for which SEDs were extracted (Fig. [2]). 



voked to mitigate this effect (Bend o et alll2010b . The beam sizes 
at 250, 350 and 500 /mi are -18.1", -24.9" and -36.6", respec- 
tively. We estimate the fluxes to be accurate within 20%. 

The MIPS 24 and 70/mi data (from MIPSGAL. ICarev et al.1 
2009 |) as well a s the I RAC 3.5 to 8 /mi (from GLIMPSE, 
Churchwell et al.l l2QQ9h and the PdBI 3.2 mm co ntin- 
uum data were first pres ented in Beuther ^t al.l (120051) and 
Beu ther & Steinackerl (120071) . On the long- wavelength side, 
we use the SCUBA 850 u rn data from the SCUBA archive 
(iDi Francesco et all 120081) and the 1.2 mm continuum flux 
measurement observed with the IRAM 30 m telescope 
(Beut heret all 120021) . The accuracy of the flux measurements 
at (sub)mm wavelength is estimated to be correct within -15% 
and for the 24 /mi within 20%. Because no sources are detected 
in the IRAC bands for most of our targets (except of the 
HMPO), the 3<x upper limits of the four Spitzer IRAC datasets 
are 0.05 mJy at 3.6 and 4.5 /mi, 0.13 mJy at 5.8 /mi and 0.15 at 
8 /mi. 

3. Results 

3.1. General structure of the complex 

Figure [T] gives an overview of the whole complex from mid- 
infrared via far-infrared to (sub)mm wavelengths. The covered 
region of approximately 4' x 9' corresponds at the given dis- 
tance to an area of ~42pc 2 with a length of ~10pc. At 70/mi, 
the linear spatial resolution is -20000AU. Globally, one can 
clearly distinguish the different evolutionary stages and the cor- 
responding spectral behavior in these images. We will discuss 
the global properties first, spectral energy distributions of exem- 
plary sources will be presented in Sect. I3.2I 

While the HMPO and IRAS source in the north of the field 
(IRAS 18223-1243) is a strong emission source at all covered 



wavelengths, this is clearly not the case for most other parts 
of the region. For example, the filamentary structure south of 
IRAS 18223-1243 is a clear (sub)mm continuum emission re- 
gion, but it shows pronounced filamentary absorption structures 
against the bright Gal actic background at 24 iim and shorter 
wavelengths (see also Beut her et al .1120051; iBeuther & SteinackeH 
2007). Nevertheless, while the (sub)mm emission peaks in the 
filament extending south to a declination of ~ -12°50'00" 
show o nly absorption in th e IRAC bands up to 8 /mi wave- 
length (Beuthe r~et all [2005). they are weak emission sources 
in the Spitzer 24 /mi band. Combining this weak mid-infrared 
emission features with outflow signatures like the "extremely 
green objects" (EGOs. lCvganowski et"aDl2008l) and CO and CS 
line wing emission favors the interpretation of these sources as 
low- to intermediate-mass protostars embedded in high-mass 
cores that will likely form ma s sive stars at the end of the 
evolution (IBeuther & Steinackerl 2007). Again different from 
these sources are the (sub)mm emission peaks found south 
of ~ -12°50'00". Almost all of these (sub)mm sources are 
either extinction features at 24 /mi wavelength or they show 
extended emission associated likely with the fore- or back- 
ground. None of these sources show clear 24 /mi emission in- 
dicative of a protostar during its formation. While some of these 
(sub)mm sources start appearing as far-infrared emission peaks 
at 70 /mi in the Herschel data, others remain in absorption or at 
least "inconspicuous" to wavelengths as long as 100 /mi. While 
more targeted studies in the past almost alway s found signs 
of sta r formation activity to ward IRDCs (e.g ., Rathb orne et al.l 
120061: IBeuther & Sridharanll2007l: iMotte et alJl2007h less biased 
Herschel studies now have the potential to identify and charac- 
terize genuine high-mass starless cores. 

From a more general point of view, these images can be con- 
sidered as templates of how young star- forming regions change 
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their appearance with wavelength: While the youngest and cold- 
est regions appear as shadows from near- to mid-infrared wave- 
lengths against the stronger diffuse background emission, back- 
ground and source emission approximately equalize at far- 
infrared wavelength (e.g., 100 /mi panel in Fig. Q]). Going to 
even longer wavelengths, the cold SEDs raise the emission of 
the cores considerably above the background, and one clearly 
identifies them as emission sources. This evolutionary picture 
strongly resemble s that of low-mass star- forming regions (e.g., 
lAndre et alJfl993b - 

3.2. Spectral energy distributions 

To better characterize the evolutionary sequence, we selected 
four sources representing different evolutionary stages: (a) the 
HMPO IRAS 18223-1243 ( Sridhara netall 120021: iBeuther et all 
2002), (b) the low- to intermediate-mass protostar embedded in 
the infrared dark filam ent IRDC 18223-3 (IBeuther et all 120051: 
IBeuther & SteinackeR 120071) . (c) a 24 /mi shadow in the south 
that becomes detectable from 70/zm onward (IRDC 18223- 
south2), and (d) a source that remains dark up to 100 /mi 
(IRDC 18223-southl). The sources are marked in Fig. [T] the 
peak positions in R.A. and Dec. (J2000.0) are (a) 18:25:10.67 
-12:42:25.8, (b) 18:25:08.61 -12:45:19.7, (c) 18:25:22.89 - 
12:54:50.7, (d) 18:25:19.15 -12:50:11.2, respectively. 

Because Herschel fully covers the peak of the SED even 
for the coldest regions, and is unique with respect to spatial 
resolution, sensitivity and robustness against saturation at the 
given wavelengths, this offers the opportunity to characterize 
this evolutionary sequence from an SED point of view in de- 
tail. Therefore we extracted the fluxes over the full wavelength 
coverage between 12 /mi and 1.3 mm wavelength covering the 
Wien and Rayleigh- Jeans parts of the spectrum. We derived 
the fluxes in all bands for uniform apertures. For IRDC 18223- 
3 and IRDC 18223-south2, the chosen aperture was based on 
the coarsest spatial resolution of the SPIRE 5 00 /mi band 
(FWHM-36.6"). For the HMPO IRAS 18223-1243, we chose a 
40" aperture to be able to combine the data with the IRAS fluxes 
at 12 and 25 /mi (Spitzer is saturated for that source at 24 /mi and 
shorter wavelengths). And for the youngest source we choose an 
even larger aperture of 60" to cover the whole clump as visible 
in Fig.[U Table [T] shows the derived fluxes and Fig. [2] presents 
the corresponding SEDs. 

The SEDs were fitted with modified Planck black-body 
functions accounting for the wavelength-dependent emissiv- 
ity of the dust. Th e assu med dust composition follows 
lOssenkopf & Henningl (Il994l) . and the assumed gas-to-dust 
mass ratio is 100. The data longward of 70 /mi can be well 
fitted with a single black-body function of cold dust and gas 
(Fig. 0. We also fitted the data for smaller apertures ignor- 
ing the less resolved 350 and 500 /mi SPIRE data points, and 
the derived temperatures agree well with our presented fits. 
Therefore, on the scales resolvable by these data for the cold dust 
and gas components, single component fits are still adequate. 
Furthermore, detections at 24 /mi and shorter wavelengths for 
IRDC 18223-3 and IRAS 18223-1243 indicate additional inner 
heating sources. This implies a temperature gradient through- 
out the inner core for these more evolved regions. However, 
a sophisticated radiative transfer modeling is out of the scope 
of this letter. Although nominally one can fit a second compo- 
nent to this warmer gas (Fig. [2]), only the luminosities of these 
components are useful parameters, while masses and temper- 
atures are poorly constrained because of rising optical depth 
at shorter wavelengths. It should be noted that in particular 



Table 1. Source fluxes in [Jy] 
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The SCUBA 850 jum data and the 1.2 mm data are taken from 
Di Franc esco et al.l d2008h and [Beuther et al. (2002), respectively. 
^RAS flux; 2 Artifact (see Fig. QJ; 3 Beuther & Steinacker (2007), 
4 Because the source is embedded in absorption filament, we selected 
smaller aperture just covering the emission source. 5 Because of ab- 
sorption close t o the source, here we a pplied PSF photometry based on 
the template by Poglit sch et al.l (2010). 6 We corrected for the mapping 
artifact visible around the source in Fig.[T] 



for IRDC 18223-3 this additional luminosity only barely con- 
tributes to the total lu minosity, consistent with pre-Herschel fits 
(IBeuther & Stei nacker, 2007) . The gas masse s calculated from 
the fits using the Ossenkopf & Henningl (1 19941) dust model with 
thin ice mantles can be considered as lower limits. Hence, we 
also calculated the gas masses from the 850 /mi data assuming 
optically thin dust continuum emission, a dust spectral index 
fi = 2 (corresponding to a dust absorption coefficient /C85o^m_~ 
0.8 cmV 1 ) resem bling the general ISM dHildebrandl 1 19831: 
Hu nter et all 120001) . and the temperatures from the cold com- 
ponent fits. Another quantity to differentiate evolutionary stages 
is the ratio of bolometric to submillimeter luminosity L/L sm that 
is usually applied to low-mas s protostars but is also discussed 
in the high-mass regime (e.g.. I Andre et aljfl993h iMueller et al.l 
2002). Table [2] lists the derived parameters. 



Table 2. Parameters for the cold (c) and warm (w) components 
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18223-3 


18 
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36 


54 


11 
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18223-s2 


16 
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157 
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21 


8 


169 


18223-sl 


18 


142 
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- 


36 


10 
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T c , M c , L c and L w are fit results. L sm is the submm luminosity (longward 
of 400 //m, and M 850 is the mass in the derived from only the 850 //m 
flux (see main text). 



4. Discussion and conclusion 

These regions represent a potential evolutionary sequence where 
IRDC 18223-southl resembles the youngest detectable stage 
with a large emission peak of cold dust continuum emission 
on the Rayleigh- Jeans tail of the SED longward of 160 /mi. 
However, at 100 /mi the source is lost in the background, and 
at shorter wavelengths it remains an absorption shadow. The 
source IRDC 18223-south2 appears a bit more evolved because 
one detects an emission source from 70 /mi to longer wave- 
lengths. Nevertheless, the region is still young and cold with- 
out any emission at 24 /mi. Continuing on the evolutionary lad- 
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Fig. 2. SEDs of the four discussed example regions shown in Fig.Q] The left panel presents an already more evolved HMPO, and 
going to the right, the regions become progressively younger with IRDC 18223-southl being only detected longward of 100/zm. 
The stars with arrows mark upper limits, whereas stars only give estimates of the existing background flux. 



der, IRDC 1 822 3-3 exhibits 24 /mi emi ssion but remains dark in 
the IRAC bands. lBeuther & Steinacked (120071) called this a high- 
mass core with an embedded low- to intermediate-mass protostar 
potentially becoming massive at the end of the evolution. Finally, 
at the northern end of the filam ent IRAS 18223-1243 comprises 
all features of a typical HMPO dSridharan et alll2002h . 

This sequence is also reflected in the SED fitting. For the 
cold sources, the single black-body fits result in low temper- 
atures between 16 and 18K, luminosities of only a few hun- 
dred L© and gas masses of several hundred M . The luminosi- 
ties of these cold sources are not related to internal heating, 
but are mainly produced by the external radiation field. The 
low temperatures are consistent with previous estimates toward 
infrared dark clouds based on spectral line observations (e.g., 
ISridharan et ail 120051: IPiilai et al] 120061) . We note that our fit- 
ted temperatures are also in the range of temperatures d erived 
from COBE for the general ISM (e.g. jReach etaDll995l) . On 
the given scales (>0.65pc) at these early evolutionary stages, 
the gas and dust clumps still have temperatures comparable to 
those of the general ISM. For the most evolved HMPO, the 
fitted cold temperature is warmer than for the other sources. 
Furthermore, the luminosity already indicates the presence of 
an HMPO, while the gas reservoir is still large, allowing for fur- 
ther accretion. While the total luminosity com pares well with 
olde r fits that were only based on IRAS data (Sridh aran et all 
2002), the cold temperature of the new fit is lower than those 
only based on IRAS data (28 versus 50 K, respectively). This is 
because IRAS alone did not sample the peak of the SED well 
and hence overestimated the temperature. Therefore sampling 
the full SED is important even for such evolved sources to de- 
rive reasonable temperature estimates. 

In the framework of the L/L sm ratio, the HMPO is above 100 
where 200 was defined as the border betwe en class and class 
I for low-mass sources (lAndre et all 1 19931) . While this cannot 
easily be translated to high-mass star-forming regions, it shows 
that even the HMPO is still in a young evolutionary stage (see 
Sridhara n et al.ll2002h . The three infrared-dark sources all have 
exceptionally low L/L sm ratios around 10, lower than the high- 
mass cores reported by Muell er et al.1 ([2002), and amon g the 
lowes t values so far reported in the literature (e.g. J Young et al.1 
2003). This supports the extreme youth of these sources. That we 
do not see a marked difference in L/L sm for the three infrared- 
dark sources indicates that - in the evolution of the massive cores 
- first signatures for star formation do not come from a global 
warm-up of the bulk of dust and gas. More evolved sources 
like the HMPO 18223-1243 in turn have multiple T components 



which show up at mid-infrared wavelengths, but most of the cold 
dust can still be reasonably well fitted by a higher but single tem- 
perature. 

While we present here one of the first sources observed with 
PACS and SPIRE, there is tremendous potential in this kind of 
observations. For the future we anticipate full radiative transfer 
modeling of such regions, which will result e.g. in more detailed 
temperature maps and density structures over the entire field. 
Furthermore, the Herschel key project EPOS (PI O. Krause) 
contains more than 40 regions in young evolutionary stages. A 
combined analysis of the whole sample should result in a robust 
characterization of the physical properties and early evolution of 
high-mass star formation. 
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